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Abstract 

In this paper we are interested in the semi-classical estimates of the 
spectrum of the Neumann Laplacian in dimension 3. This work aims 
to present a complementary case to the one presented in the paper 
of Helffer and Morame in the case of constant magnetic field. More 
precisely, in the case when the magnetic field is variable and under 
the most generic condition for which boundary localizations can be 
observed, we prove a three terms upper bound for the lowest eigenvalue 
and establish some semi-classical behaviour of the spectrum. 

1 Introduction and main results 

Let Q be an open bounded subset of M 3 with smooth boundary. For a vector 
potential A e C°°(n,M 3 ) and h > 0, we introduce 

(3 = V x A 

and the quadratic form defined for all ip 6 H C) by : 

QM) = [ \(ihV + A)ij\ 2 dx. 
Jn 

One of our interests will be the lowest eigenvalue denoted by Aj'(A) of the 
associated self-adjoint operator i.e. the Neumann realization of (ihV + A) 2 
on Q denoted by P£. The question of analyzing the behaviour of Aj(A) as 
h tends to appears in the theory of superconductivity when studying the 
third critical field of the Ginzburg-Landau functional (see |FH09| ) and also 
in the theory of phase transition of liquid crystals (see |HP07| ). In this last 
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topics, it turns out that the case when the magnetic field admits constant 
strength really plays a particular role. In this paper, we will treat the case 
when both direction and strength of the magnetic field may vary, what can 
happen in the superconductivity case. In fact, until now, only the case with 
constant magnetic field was studied in details. Before analyzing the general 
problem, we need to recall some properties of a simplified model on the half 
space which were first established in |LP00| and then in |HM02[ |FH 09|. 

Model in M. 3 + 

Let us consider the case of the constant magnetic field in 

= {( r)S) t) el 3 :t>0}. 

The angle between the magnetic field (3 and the plane t = is denoted by 9\ 
and so, up to a rotation, we can assume (in the case when the norm of (3 is 

1): 

(3 = (0, cos 9, sin 9), 

with 

An associated vector potential is given by : 

A = (tcos9 - ssin#,0,0) 
and we consider the Neumann realization on of 

(1.1) H{0) = (D r + V e ) 2 + D 2 + D 2 , 
where 

(1.2) V e = tcos9 - ssin#. 

The bottom of the spectrum of this operator is denoted by : 

(1.3) mf a (H(9)) = a{9). 

Atlthough we work in 3D, we will have to consider the Neumann realization 
on denoted by H{6) with domain : 

D(H(9)) = {tpe L 2 (M. 2 + ) : H(6)tp G L 2 (M. 2 + ) and D t tp = on t = 0} 

and defined by : 

(1.4) H(9) = D 2 + D 2 + V 2 . 

It has been proved in |LP00| (see |FH09| and the references therein), that : 
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1. inf a{H{9)) = o{6) for 6 e ] 0, § [, 

2. a is analytic on ]0, |[, strictly increasing on [0, |] and satisfies : 
< <t(0) < 1 and a (f ) = 1. 

3. For 6 e ]0,f [ : 

(1.5) inf<7 esfl (tf(0)) = 1. 

4. cr(0) is a simple eigenvalue of -ff"(0) associated to some positive and 
L 2 -normalized function uq (see also [RS78]). 

In particular, the minimal energy of 1~L{6) is obtained when the magnetic 
field is tangent (0 = 0). Using a rescaling, the bottom of the spectrum when 
the field admits \\(3\\ for module is given by ||/3||cr(0). 
For x G dfl, we let : 

(1.6) $(z) = *(p(x))\\/3(x)l 

where 9(x) is defined by : 

\\P(x)\\sm9{x) =(3-v(x) 
with v{x) the inward pointing normal at x. 

Lu-Pan's result : main term of the asymptotics 

We recall the asymptotics which was first established by Lu-Pan (cf. [LP001 
IFH09] ) : 

Theorem 1.1 The following asymptotic expansion holds : 

X'l(A) = min(inf ||/3(x)||, inf ${x))h + 0(/i 5/4 ). 

Remark 1.2. 

1. In fact, until now, the proof of the upper bound was just sketched (see 
|LP00l Theorem 5.1]). We will see that, if some generic conditions are 
satisfied, the remainder for the upper bound is better than 0(h 5 ^). 

2. The Stokes formula implies that : 

(3-uda= / V -(3dx = 0, 
an Jn 

and so there exists a point where the magnetic field is tangent to 
the boundary. Consequently, when the magnetic field has constant 
strength, the first term in the asymptotics is 0q||/3(xo) \\h. 
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3. The case when the magnetic field is constant has already been studied 
in details under generic assumptions in |HM04| where a two terms 
asymptotics is proved with a second term of order /i 4 / 3 . 

4. When the magnetic field just admits a constant strength, one or two 
term asymptotics are established in the papers |HP07i IHP081 Ray09b| . 



5. Finally, let us mention that this kind of asymptotics has already ap- 
peared in dimension 2 in the constant case (see |BS98l IBPT981 ILP991 
IdPFSOOl IHMOll rFH06| ) and that we have obtained a two terms expan- 
sion in |Ray09a| when the magnetic field is not uniform. 



Main assumptions and ideas 

Let us state our hypotheses on the magnetic field in this paper. Our estimates 
are only interesting when the following condition is satisfied (cf. Theorem 
tDJ): 

(1.7) inf/3<inf||/3||. 

ail i I 

This condition is also called the surface superconductivity condition (cf. 
|LP00l |FH09| ). It implies a localization of the groundstates on the boundary 
when h — > 0. Hence, on the boundary, we wish to deal with some generic case 
in the same way as we have done in 2D (see |Ray09a| ) ; so we will assume 
first that : 

(1.8) (3 admits a minimum at xq 

and that at this point the field is neither tangent, neither perpendicular to 
the boundary, i.e : 

(1.9) < 6(x ) = 0i<\- 

Our point is to obtain the semi-classical behaviour of eigenvalues in the bot- 
tom of the spectrum of In order to do that we will construct an operator 
whom spectrum is close to the one of P£. From this, it will follow a quasimode 
construction associated to the approximated operator and simultaneously we 
will get, in some sense, the candidate to a three terms asymptotics (in powers 
of h) of the lowest eigenvalues. To define the different terms in the asymp- 
totics, each one corresponding to a step of the quasimode construction, we 
need to introduce a few invariants mixing the geometry of the domain, the 
properties of the magnetic field and spectral quantities attached to H{6). We 
first define the direct orthonormal basis ( Iq , Tq , Vq) attached to the point x : 
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1. Vq is the inward pointing normal at Xq 

_> (3(x ) - {(3{x ) ■ 

||/3(x ) - (/3(x ) • I'oHII 

3. Z = tq x ^o- 
Let us now present our main results. 

Second term of the asymptotics 

Our first invariant depends on the the second fundamental form at x denoted 
hyK: 

K n = K(X, X), K 22 = Kffi, K 12 = K(X, t£). 
It also depends on the first derivatives of (3 at xo : 

V = V < \\ a f mp To > -"0 - K n cos(0), 
POo)|| 

with 9 = 9\. At x we define the quantity C^' K (xo) by : 

2/3 < V e (st - —s 2 )u e , u e > +rj < t 2 V e u e , u e > +^— - / t\D s u g \ 2 + < H K u e , u B >, 
2 sin 9 J t>0 

(1.10) 

where is defined in (11.21 ) and 

(i.n) r = r W = -§|, 

with 

(1.12) C = C(0) = cos 9(7(6) - sin 0a' (0) , S = S(9) = sin 9 a (6) + cos 0a'(0) 
and 

(1.13) H K = (K n + K 22 )d t + 2tK u V e 2 + 2tK 22 D 2 + 2tK 12 (V e D s + £> a V&). 
Remark 1.3. 
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1. We can notice that S(9) > 0. Indeed, we know (see |FH09| ) that 
cr'{9) > 0, a (9) > and we have assumed that 9 ^ 0. 

2. An integration by parts provides that 

< t(V e D s + D 8 Vo)u , u e >= 
and so C^ ,k {xq) is real. 

3. We will see that the construction of C^ ,k {xq) as part of the second 
term of the asymptotics uses that x$ is a critical point of (3. 

■ 

We can now state a theorem dealing with the second term of the asymptotics : 

Theorem 1.4 Under Assumptions ( li.gj) and $1.9]) . there exists D > and 
ho > such that, for < h < h , there exists at least one eigenvalue \i of 
P\ such that : 

|/i - 0(x o )h + C^(x )||/3(a;o)ir /2 / i 3/2 )| < Dh\ 
where C^' k {xq) is defined in $1.10]) . 

Then, to continue the asymptotic expansion, we assume the following generic 
condition : 

(1.14) /3 admits a non-degenerate minimum at Xq. 

Third term of the asymptotics 

We denote by &p the Hessian matrix at x , namely : 

6/3(r, a) = \d 2 r P{x Q )r 2 + ^J(x )(sr + rs) + \d 2 J(x )s 2 

We can now introduce a fundamental operator linked with the behaviour of 
the third term of the asymptotics : 

(1.15) 6^ = 6^,^). 

sin 9 

By the non-degeneracy Assumption (11.141) . this operator is an harmonic os- 
cillator whose spectrum is explicitly known. We can now state the main 
result of this paper : 
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Theorem 1.5 Under Assumptions Ijl.l4\ ) and U.9\) . there exists d £ R such 
that for all n£N*, there exists D n > and h n > such that, for < h < h n , 
there exists at least one eigenvalue fi n of P\ such that : 

\li n - 0(x o )h + C^(x )||/3(x )ir/ 2 /i 3 / 2 + ( ln (&?) + d)h 2 )\ < D n h 5 / 2 , 
where 7„(@ / g) is the n-th eigenvalue of 
This theorem admits the following easy corollary : 

Corollary 1.6 Under Assumptions jl.l4\ ) and (COjj, there exists d G E, 
Di > and hi > such that, for < h < hi : 

X h i(A) < P(x )h + C^ K {x Q )\\f3{x )\\ l ' 2 h^ 2 + ( 7l (e^) + d)h 2 + Dih 5 / 2 . 

An interesting point would be to get the corresponding lower bound. 



Organization of the paper 

In Section 2, we study some spectral properties of H(9). In Section 3, we 
consider a change of coordinates and see how the magnetic field and the 
metrics are transformed. In Section 4, we explain how we get, using Tay- 
lor approximations, a simplified operator. Finally, in Section 5, we use the 
simplified model to construct a quasimode and prove Theorem 11.41 and 11.51 



2 The family of operators H(6) 

This section is devoted to the study of a family of operators H(6) which 
plays a fundamental role. Let us explain why we are reduced to this family 
of operators in . Considering the operator ( 11.1ft and performing a Fourier 
transform in the variable r, we are led to the r-dependent family of Neumann 
realizations : 

H(9,r) = (r + V e ) 2 + D 2 + D 2 . 

Then, under Assumption fll .OR and making a translation in the variable s, 
H(9,t) is unitarily equivalent to H(9,0) = H{6). 

In the next subsection, we study the exponential decay properties of the 
solutions of equations in the form : 

h(9)v = /, 

where : 

(2.16) h{9)=H{9)-a{9). 

This kind of equations will indeed appear many times in the final quasimode 
construction. 
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2.1 Exponential decay 



We begin to notice that, in view of (11.51) and using Persson's Theorem (cf. 
|Agm82[ IFH09j ). when < 9 < f and for any < a < y/l - a (9), we get : 



Let us denote 

L 2 exp (M 2 + ) = if e ■ 3a > 0, e <*(*+<*»/ G L 2 (R 2 + )} 

and 

= {/ G = 3a > 0, e a < t+< '»/ G tf 1 ^)}, 

where < s >= (s 2 + 1) 1/2 . 
Thus, we have : 

(2-17) u e G ^ P (M 2 + ). 

Let us now state the main proposition of this subsection. 

Proposition 2.1 Let f G Ll xp (M.+) such that < f , ug >= 0. Then, if v 
denotes a solution of 

(2.18) h(9)v = f, 

then 

v G Hl xp (R 2 + ). 

Proof. 

We begin to prove the control in the variable s. Let us introduce a smooth 
cutoff function x satisfying 

X(s) = 1 if \s\ < 1 and x(s) = if \s\ > 2. 

For n £ N*, we consider the (bounded and inversible) multiplication oper- 
ator by e~ ex "( s ) <s> where Xn(s) = x(ra _1 s). Let us verify that it preserves 
D(H(9)). 

For rj) G D{H{6)), we have : 

H{9){e- €X " is)<s> ifj) = e~ eXn{s)<s> H{9)i)-ex' n {s) < s > D s ^-ex n {s)—^—D s ij. 

As D s ip G L 2 (R 2 + ), we deduce that H(9)(e- eXn ^ <s> i)) G L 2 (R 2 + ). Moreover, 
e ~fx»( s )< s >^ satisfies the Neumann condition and so : 

e - eX n( s )<s>^ e d{H{9)). 
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We now use a Grushin method and we introduce fj(9) defined on D(H(9)) xC 
by: 

h(9) u e 
< ■ , u e > 



HO) 



It is obvious that this operator is inversible. (12 . 18H is equivalent to : 



V 




"/" 











if we take for v the unique solution orthogonal to uq (all the other solutions 
will satisfy the decay properties because uq does). Letting v = e eXn ^ s ' <s> v 
and / = e eXn ^ <s> f , we can rewrite : 



V 




" / " 







_ 



where 



g£XnW<S> Q 

1 



e -£Xn(s)<s> Q 

1 



We can notice that D(S)(9) e,n ) = D (f)(9)). An easy computation provides : 

W'™ = [ < . ; e - e ^ <s> u e > 

We can write, using the exponential decay of uq (see (12.171) ). that : 

f)(9) e ' n = f)(9) + eF n , 

where F n is bounded relatively to f}(9) uniformly with respect to n. The 
conclusion is then standard ; for e small enough f)(9) e ' n is inversible and 
there exists C(e) > such that for all n e N* : 

\\(f3(9Y' n r 1 \\<C(e). 

This leads to the inequality : 

INIl2(r2 ) < C(e 



Thus, we deduce : 



L 2 (R2 )• 



,fXn(«)<S> 



||e^( s ) <s> i;|| L2(R2+) <C(e)||e e j \\ L ^ +) . 

The dominated convergence Theorem shows that the r.h.s converges when 
n — ► +oo to || e t<s> f ||z,2( R 2). So, we infer that e eXn ^ <s> v weakly converges 
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(after subsequence extraction) in L 2 (R+) to some function h. The conver- 
gence in the sense of distributions and again the dominated convergence 
theorem proves that : e e<s> v = h E L 2 (R+). Consequently, the decay with 
respect to s is proved. 

We now explain the decay with respect to the variable t. We proceed exactly 
in the same way, but we have to deal with the Neumann condition. For 
n E N* , we consider a smooth cutoff function rj n such that : 

f if < t < \ 
rj n (t) = < 1 if 1 < t < n . 
{ if t > 2n 

We observe that the (bounded and inversible) multiplication operator by 
e ^„(t)t preserves D(H{6)) ; then, it remains to use the same analyis as previ- 
ously by using this time the decay of ug with respect to t and the proposition 
is proved after having noticed (by integration by parts) that : 

q e (e a/2{t+<s>) v) < +oo, 

where qe is the quadratic form associated to H(9). 

m 

Let us denote : 

H? xp = {fe L 2 (R 2 + ) : V(*, k) E N 4 DiD k t f E L 2 exp (Rl)}. 

It is clear that H^. p C S (R + x R) . Proposition 12.11 permits to prove the 
following proposition : 

Proposition 2.2 The groundstate uq of H(Q) belongs to . 
Proof. 

We first recall that ug E „(R+). Then, let us take the derivative with 
respect to s : 

h(6)d s u e = 2 sin 9V e ug E L 2 exp (R 2 + ). 

Proposition 12.11 provides D s u e E Hl xp (R? + ). A recursion gives that, for all 
m E N : 

D?u e E H l exp (R 2 + ). 
Then, we come back to the equation : 

D 2 t u e = (-D 2 - Vg 2 + a(6))ug E L 2 (R 2 + ). 



10 



Taking the derivative with respect to s, we find : 



D 2 t DTu e G L 2 (R 2 + ) 



and differentiating with respect to t, we deduce : 




We can now state the other important corollary : 
Corollary 2.3 Let g G Hf xp and f G D(H(9)) such that : 

h(9)f = g. 



Then, f G H™ p . 
Proof. 

The proof is essentially the same as the one of the previous corollary. 

■ 

2.2 A few properties of uq 

We begin to establish some properties of ug related to the tangential coordi- 
nate s. 

2.2.1 A momentum formula 

The definition of ug provides : 



We differentiate with respect to s this identity and obtain the formula : 



Therefore, taking the scalar product with ug, we deduce the lemma : 
Lemma 2.4 

For all g]0, |[, we have : 



H(6)u e = a(9)u e . 



(2.19) 



h(9)d s ug = 2sin9VgUg. 
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We again take the derivative in ((2. 191) with respect to s to find : 

(2.20) h(9)d 2 u e = 2sme{2V e d s u 6 - sin6u e ) = 2sin9(d s V e + V e d B )u g . 
If we take the derivative two more times, we get : 

h(9)dtu e = 12 sin 2 6d 2 ug + 8 sin 9Vgd 3 s u e . 
We deduce the lemma : 
Lemma 2.5 We have : 

Vgd 3 ugu e dsdt = ^ sin 9 J d^UgUgdsdt. 

2.2.2 Another formulas 

We will meet the question of finding g s.t : 

(2.21) h(9)g = -td s u e . 
We let : 

(2.22) L s = -d 2 s + sin 2 9 s 2 . 

An easy computation gives the following commutator formula : 
L S H(9) = H(9)L S + 4 cos 9 sin 9td s . 

We infer : 

h(9)(L s ug) = — 4 cos 9 sin 9td s ug. 
Hence, we can define a solution of ( 12.211) by : 

(2.23) f = (2 sin(2^))^ 1 (sin 2 9s 2 u e - d 2 8 u g ). 
Finally, we have : 

L s h(9)f = h(9)(LJ )+4cos6sin9td s f . 
We take the scalar product to get : 

< L s h(9)f , ug >=< 4cos9sin9td s f , u e > . 

But, we have : 

< L s h(9)fo, ug >= — < L s {td s ug) , ug >= — sin 2 9 J ts 2 d s ugug dsdt 



= sin 9 J tsu d dsdt. 
So, we deduce the following lemma : 
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Lemma 2.6 We have the identity : 

< 4cos9td s fo,uo >= sin/9 J tsu 2 e dsdt. 

Finally, differentiating twice the equation (12.211) satisfied by fo, we obtain : 

Lemma 2.7 We have, on one hand : 

h(d)d s f = -td 2 s u e +2wx9Vef , 

< -td 2 u e + 2 sin 6V e f , u e >= 
and on the other hand : 

h(6)d 2 J = -tdlue + 2 S m9(- S m9f + 2V e d s f ), 

< -tdlue + 2 sin 0(- sin 6f + 2V e d s fo),u e >= 0. 

2.2.3 Properties of H(9, p) 

The main ingredient in the following analysis is the so-called Feynman- 
Hellmann formula. We introduce a parameter p which will permit to play 
with a scale invariance. We define : 

H(9,p)= 1 -(D 2 + D 2 ) + pV 2 
P 

and consider : 

(2.24) H{e,p)u 0tP = <T(6)u e , p , 

where 

u e , p (s,t) = u e (p~ l/2 s,p~ 1/2 t). 
We can notice that H(9, 1) = H{9). 

2.2.4 Feynman-Hellmann with respect to p 

We take the derivative with respect to p : 

(H(9,p)-a(9))d p u = -d p Hu, 

and 

d p H = -p- 2 (D 2 + A 2 ) + p' l / 2 {r + y/pV e )Ve. 
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We use fl2~24l to obtain : 

(2.25) (HiO, p) - a)d p u = ^ - 2V 2 u 9>p . 

Multiplying this last equation by uq. p and integrating, one recovers the 
called Virial Theorem : 

(2.26) [ (\D s u e \ 2 + \D t u e \ 2 ) dsdt = [ \V e u e \ 2 dsdt = ^ . 

Jt>0 Jt>0 * 

Remark 2.8. 

If we only perform a rescaling in the variable s, we find the identity : 

(2.27) J \D s u e \ 2 dsdt = -sin6 J sV e u 2 9 dsdt. 



2.2.5 Feynman-Hellmann with respect to e 

We take the derivative of (12.241) with respect to 6 : 

(H(6,p) - a{9))d e u = a'(6)u - d e Hu, 

and 

d e H = 2pV e d e V e . 

We find : 

(2.28) (H - a)d e u = a'(6)u e - 2pV e d e V e u e . 
Taking p — 1 and multiplying by Ug, we get : 

(2.29) 2/ V e ^\u e \ 2 dsdt = ( r'(9), 

Jt>o 

which also reads : 

(2.30) a'(6):=-2 (cos9t - sin#s)(sin#t + cos 6 s)\ue(s,t)\ 2 dsdt . 

Jt>o 

Remark 2.9. 

Combining (12T291 and (12T26D . it follows : 

(2.31) J tV u 2 e dsdt = and J sV e u 2 e dsdt = ~^y^ 
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2.2.6 Consequences 

A first identity 

We get the following lemma : 

Lemma 2.10 Defining 

(2.32) h{0)= [ (t- sT(9))V e \u e \ 2 dsdt , 

Jt>o 

where T(6) is defined in U.ll\) . we have : h(0) = 0. 
Proof. 

One multiplies (12T251 ) by a' and (I2T28!) by - p to obtain : 

(2.33) (H(9, p) - a{9))w B , p = -2a% 2 u e , p + 2aV g dgV g ug tP , 
where 

(2.34) w e , p = cr'd p ug tP - ^-od e u e ,p. 

We take p = 1, multiply by ug, integrate and Lemma 12.101 is proved. 

■ 

We will denote : 

(2.35) w = a'dpiig^ - ad u 0>1 . 

A second identity 

We now apply the operator 

(2.36) j(e )P )=<T , d p -~d e 

to (12.331) . take p = 1 and we obtain : 

j(6, p)(H(d, p))w B , p + (H(0, p) - a)j(6, p)wg, p 
= (-2a'V e 2 + 2aV e d e Ve)we,p + j{9, p)(-2a'V 9 2 + 2aV e dgVg)u e>p . 

Then, on the one hand, we get : 

(2.37) j(6, p)(H(0, p)) = -a\D 2 + D 2 t ) + a'V 2 - 2aVgd e V 6 , 

and one the other hand : 

j(6,p)(-2a'V 2 + 2aV e dgV e ) = -adg{-2a'V 2 + 2aVgd e V e ) 
= 2a"oV e 2 + 2<ja'Ved e Ve - 2a 2 (d e V e ) 2 - 2a 2 V e d 2 e V e . 

Taking the scalar product with uq, and using (12.291) and (12.261) . we deduce 
the lemma : 
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Lemma 2.11 We have the identity : 

(2a% 2 - 2aV d e V e )w o u dsdt + / (ad e V e - a'V ) 2 u 2 e dsdt 

Jt>o 



a 2 a" + a 3 
2 ' 



A third identity 

If we take the derivative of (12.331) two times with respect to s , we get : 

h(9)d s w = 2 sin 9V e w - 2V e u e + 2(Cs + St) sin 0u e - 2(Cs + St)V e d s u e 
and then : 

h(9)d 2 w = 4 sin 6V e d s w - 2 sin 2 8w - 2(Cs + St)V e d 2 u e 
+ AC sin 9u e - 4CV 9 d s u e + 4 sin 6(Cs + St)d s u e 

and we deduce the lemma : 

Lemma 2.12 We have on one hand : 

2sin6> < Vowo, u e > +2sin^ < (Cs + St)u e , u e > -2 < (Cs + St)V e d s u e , u e >= 

and on the other hand : 

< 4 sin 6Ved s wo — 2 sin 6w - 2(Cs + St)V e d 2 u e , u e > 
= -sin9< (2C + 4(Cs + St)d s )u e , u e >= 0. 

Two more identities satisfied by / and w 

If we rescale the equation (I2.2ip satisfied by / , we get : 

(H(9,p)-o-(0))f o , p = -td s u e , p , 

where 

f 0iP (s,t) = f (p- 1 / 2 s,p-^ 2 t). 
We now apply the operator j, take p = 1 and use the identity ( 12.371) to get : 

h(e)j(f , P ) + a'td s u e + {2a'V 2 - 2aV e d e Vg)f = -td s w . 

Taking the scalar product with u$, we find : 

(2a'V 9 2 — 2aVedgVe) fougdsdt + J td s woUgdsdt = 0. 

Moreover, noticing that 

— j (^J td s u 2 ep dsdt^j = 2 J tug iP wg iP dsdt, 

we get the lemma : 

16 



Lemma 2.13 We have the following cancellations 



J (2a'V e 2 - 2aV e d e V e )f u d dsdt = 

and 

J td s WQUgdsdt = 

Remark 2.14. 

It follows from Corollary 12.31 that all the solutions of equations of the form 
h(6)f = g that we have met in this section belong to H™. 



3 A choice of coordinates near the boundary 

This section deals with a system of coordinates near the boundary. Indeed, 
it will be useful to be reduced to a simplified model on an half space. 

3.1 A general choice of coordinates 

Let us assume that £ dVL. In some neighborhood V of 0, we take the 
coordinates (yi, y 2 ) on dQ (via a map C 3 denoted (j) and which will be defined 
precisely in the following). We denote 2/2)) the inward pointing 

normal at the point 4>~ 1 (yi,y 2 ) and we define local coordinates in V : 

$(j/i,?fe,2fe) = 2/2) + y3"((i>~ 1 (yi,y2))- 

More precisely, for a point x £ V, 2/2) is the projection of x on dQn V 

and V3 = t = d(x, dQ). 

Transformation of the magnetic field 

We now want to determine the new vector potential and magnetic field in 
these new coordinates. Let us introduce the 1-form uj : 

uj = A\dx\ + A 2 dx 2 + A 3 dx 3 . 

We have : 

u = Axdyi + A 2 dy 2 + A 3 dy 3 . 
We can express dxi as a function of (dyj) and we get : 

A = D y $- 1 (A($- 1 (y))). 
17 



Then, in order to find the new field, we write : 

duj = (V x A)idx 2 A dx 3 + (V x A) 2 dx\ A dx 3 + (V x A) 3 dxi A dx 2 . 
The comatrix formula provides : 
(3.38) /3 = det(D$)- 1 ((D$))'/3, 

where 

P = V y x A. 

Metrics in the new coordinates 

The euclidean metrics is 

go = dx\ + dx\ + dx\. 
In the new coordinates (2/1, ?/2j 2/3) , we have : 

9o = ^gijdyk^dyj, 

where (gkj) is the matrix We denote (<7 fcj ) the inverse matrix 

of (gkj)- 

If the support of u is sufficiently concentrated near Xq, we have : 
(3.39) 

Qa( u ) = / + A)w| 2 dx = / \(ihV y + A)M|^ i(D$) |detD$- 1 |(iy 1 

where 

I] 9 kj (ihV yk + A k )u (ihV yj + Aj)u. 

We will denote by Pj£ the associated Neumann realization on L 2 (| det D$ _1 |dy), 
namely : 

(3.40) Pi = \ <let D^\^(ihW j + A^det D^~ l \g kj (ihW k + A k ). 
3.2 Coordinates (r, s, £) 

We now introduce normal coordinates near Xq. (tq, 1$) is an orthonormal basis 
of the tangent plane at x$ ; we denote by (r, s) the corresponding coordinates. 
Then, it is standard, using the exponential map near xq, that (r, s) define a 
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local parametrization of the boundary and they are called normal coordinates 
(cf. |Laf96j). So, there exists an open set S of M 2 and a diffeomorphism (f> 
such that : 

<f>:W X0 ^S,<f>{x) = (r,s). 
Moreover, if G denotes the first fundamental form of dfl, we have : 

G = / 3 + Gi(r, s) + 0(r 3 + s 3 ), 

where G\ is a quadratic form in r and s. 

Letting t(x) = d(x,dQ), we define the system of coordinates on the bound- 
ary : (2/1,^/2,2/3) = (r,s,t). The metrics go can be expressed as : 

g = dt®dt + G- 2tK + t 2 L, 

where K, L are respectively the second and third fundamental forms on dQ. 

Taylor expansion of the metrics 

We let 

' ATn K 12 



K, 







K\2 K 2 2 

We write the followimg Taylor expansion : 

K = K + K 1 (r, s) + 0(r 2 + s 2 ] 





u 12 





u 12 


U22 















where K\ is linear in r and s. 

In general, for a 2 x 2 matrix U, we will denote U° the matrix 



U° 



Thus, we can write : 

g = h- 2tK° + G?(r, s) - 2tK% + t 2 L° + 0(r 3 + s 3 ) 
= I 3 - 2tK° + G°(r, s) + R(r, s, t) + 0(r 3 + s 3 + t 3 ) 
= ^a PP + 0(r 3 + s 3 + t 3 ), 

where R is homogeneous of degree 2 and of partial degree 1 with respect to 
r and s. 
Remark 3.1. 

In the following, we take as convention to always denote by R the homo- 
geneous polynoms of degree 2 and of partial degree 1 with respect to r and 
s. 
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We deduce : 



and 




It follows : 



g\ 1/2 = l-tK^ + Ti(^) + R{r,s,t) + 0{r 3 + s 3 + t 3 ) 



m app + 0(r 3 + s 3 + t 3 ), 



where K, 



Kn + K 2 2 is the mean curvature and : 



m, 




1/2 



app 



The dual metrics satisfies : 



(3.41) 

g° = h + 2tK° - G?(r, s) + R(r, s, t) + 0(r 3 + s 3 + 1 3 ) = g° app + 0(r 3 + s 3 + 1 3 ). 



Remark 3.2. 

In the following, it will be convenient to denote by (resp. Gf ■ , the 
coefficient of r 2 in the coefficient with index of G\ (resp. the coefficient 
of s 2 , the coefficient of rs). 



4 An approximated operator on 

In this section, we use the local coordinates introduced in the previous sec- 
tion to construct a model operator near Xq. The one term asymptotics was 
obtained by considering the model with constant field in 1R 3 ; in the follow- 
ing, we keep more terms in the approximation of P£ to get a three terms 
asymptotics. 



In this subsection, we approximate first the metrics and then the magnetic 
field in order to get a model operator. 



4.1 The model 
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Approximation of the metrics 

With (13.401) and f!3.4ip . we are reduced to consider the operator on L 2 (m app drdsdt) , 
with Dirichlet condition on t = t , \r\ = r , \s\ = So (for to, r , So positive 
and small enough) and Neumann condition on t = : 

H app =h 2 (m app )- 1 D t m app D t 

+ (1 + 2tK n - G n - R n )(hD r + A r f + (1 + 2tK 22 - G 22 - R 22 )(hD s + A s f 
+ (2K 12 t - G 12 - R n ){hD r + A r )(hD s + A.) 
+ (2K 12 t - G 12 - R 12 )(hD s + A s )(hD r + A r ) 
+ h 2 Rx{r, s, t)D r + h 2 R 2 (r, s, t)D s , 

where the Ri(r, s, t) are linear functions in r, s and t coming from a commu- 
tator between the metrics and the derivatives D r and D s . 
This first approximation of the operator Pj£ satisfies, for if) sufficiently con- 
centrated near Xq (an explicit choice will be made in the following) : 

(4.42) 

||P^ _ u-vv^w < Ch 2 ||(|r| 3 + |s| 3 + t 3 )\D k Dji;\ || 

+ Ch\\(\r\ 3 + | s | 3 + t 3 )|^||| + C7i^||(r 4 + S 4 + t 4 )|Z^||| +C||(|r| 5 + | S | 5 + t 5 )|^|||. 

j 

Omitting the terms of order 3, we get : 
(4.43) 

H MOD = h 2 n app D t m app D t 

+ (1 + 2tK u - G n - R u )(hD r + A T ) 2 + (1 + 2tK 22 - G 22 - R 22 )(hD s + A s ) 2 
+ (2K 12 t - G 12 - R 12 )(hD r + A r )(hD s + A 8 ) 
+ (2K 12 t - G 12 - Rn)(hD s + A s )(hD r + A r ) 
+ h 2 Ri{r,8,t)D r + h 2 R 2 (r,s,t)D s . 

with 

n app = l + tK* 1 -Tr(-±)+R. 

Taylor approximation of the magnetic field 

The Taylor formula applied to (3 at x gives (at the order 1): 

P = A, + rd r (3(x ) + sd s f3(x ) + td t /3(x ) 
= (3^ + 0(r 2 + s 2 + t 2 ), 
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and with (13.381) we are led to : 

/3 = (3 + t(-K^/3 + A'°/3 + dt0(x o )) + rd r /3(x ) + sd s /3(x Q ) + 0(r 2 + s 2 + £ 2 ). 
= ^+0(r 2 + s 2 + t 2 ). 

We can notice that, on t = 0, we have = /J^. 
Let us write : 

Pr (1) = S r - e s - f *, 

j3 s ^ = cos Q x + a r + /3 s + r^t, 

Pt m = sin 6>i + 7 s + (or + A*o*, 



and 



(3 r ^ = —5r — es — ££, 



0: 



o (1) 



3(1) 



cos 6*i + ar + /3s + 
sin + 7s + Qr + /it. 



Then, we have the relations : 

(4.44) S = 5, e = e, a = a, /3 = p, To = 7> Co = C, 

(4.45) £ = £o — -^12 cos &,V = Vo ~ K\\ cos A* = A*o — -^o^ si n 

We now write the formula at the order 2 for j3 and but only when t = ; 
indeed, this will be enough in the following. So, the approximation at the 
order 2 of (5 and /? on t = can be written as : 

j3 M = -Sr -es- C~ Q rs - F^s 2 - A> 2 , 

/3 S ( 2 ) = cos 9x + ar + (3s + C rs + F s 2 + D r 2 , 

(3 t {2) = sin Q x + 7 s + (r - 2B rs - 3H s 2 - E r 2 . 



/V =-5r-es- Crs - Fs 2 - Dr 2 , 

(3 S {2) = cos 9 1 + ar + ps + Crs + Fs 2 + Dr 2 , 

(3 t (2) = sin 0! + 7 s + (r - 2Brs - 3Hs 2 - Er 2 . 
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Then, we have (cf. 13.381) . on t = : 

^(l +Tr ,|)) (/ 3 -f)/3 

and this provides a relation between the capital letters of $ and fi. Identifying 
the terms in s 2 , we get (cf. Remark l3~2l) : 

(4.46) F = T + -coaeG^, 

F = F + ^cos9G( 1 , 

H = H Q -\sm9{G( 1 + G s * 2 ). 
o 

Then, we obtain for the term in r 2 : 

(4.47) £ = A) + -cos0G£, 

D = D + - cos 9G£, 

E = E -^sm6(G\l + G£). 

and for the term in rs : 

(4.48) C = C^+- cos 9G r 12 , 

1 

11 5 



C = C + =rcos9G r 



2 

£ = 5 o -Tsin0(G™ + G2). 
Then, we choose a gauge such that : 

A t = 

and such that the Taylor approximation of A at the order 3 denoted by 
is defined by : 

i r (3) = V + P r>2 + P r , 3 

i, ,3) = o, 



23 



where 

p r 2 = (at - (s)r + (3st - ^s 2 + |t 2 
(4.49) 

p r3 = (At 2 + Bs 2 + Cst)r + (Dt + Es)r 2 + (Fts 2 + Gt 2 s + Hs 3 + It 3 ) 
P s ,2 = 8rt + est + |t 2 

p s 3 = At 2 r + Cstr + Dtr 2 + Fts 2 + Gt 2 s + Jt 3 . 
The model operator on L 2 (m app drdsdt) is : 
(4.50) 

U M = h 2 n app D t m app D t 

+ (1 + 2tK n - G n - R n )(hD r + i r (3) ) 2 + (1 + 2tF 22 - G 22 - R 22 ){hD s + i s (3) ) 2 

+ (2K 12 t - G 12 - R 12 )(hD r + A-^)^ + A s (3) ) 

+ (2ifi 2 t - G 12 - R 12 ){hD s + i s (3) )(M) r + i r (3) ) 
+ h 2 R 1 (r,s,t)D r + h 2 R 2 (r,s,t)D s 

and it satisfies : 
(4.51) 

\\H app ip - H m 4j\\ < ChJ2\\(r 4 + s 4 + t^D^W + C\\(\r\ 5 + |s| 5 + t 5 )\if;\\\. 

3 

In conclusion, we can write : 
(4.52) 

\\P& - w^vil < ch 2 IKM 3 + l^l 3 + * 3 )|£WIII 

+ C/i||(|r| 3 + |s| 3 + t 3 )|^||| + ChJ2 \\(r 4 + s 4 + t 4 )|D^||| + C||(|r| 5 + |s| 5 + t 5 

j 

4.2 Non-degenerate minimum of 

In this subsection, we express the condition of non-degenerate minimum of 
(3. Firstly, we write that x is a critical point of (3 and secondly we give 
the expression of the Hessian matrix of (3. Moreover, we will assume that 

!l/3(zo)|| = i. 
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4.2.1 Critical point of (3 

By definition of 9, we have : 

(4.53) ||/3|| sin = A- 

Let us compute the partial derivatives of 9. We have : 

d r& = ^Jpf s ™9+\\f3\\cos9d r 9 
d s p t = ^^-sm9+\\f3\\cos9d s 9 



In r = s = t = 0, we get 

c 

7 

and so : 



3 get : 

= (a cos 9\ + C sin #i) sin #i + cos 
= (/3 cos 6*i + 7sin^!) sin^ + cos9id s 9 



d r 9 = C cos 9\ — a sin 6^ 
<9 s /9 = 7 cos 9\ — (3 sin 6*i . 



Then we differentiate /3 with respect to r and s : 



At r = s = t = 0, we have : d r (3 = and d s f3 = 0. This leads to 

(a cos 9 1 + C sin 0i)o-(0i) + a'{9 x )d r 9 = 
(/3cos0i +7sin# 1 ) C r(# 1 ) + </(#i)a s 6> = 



Thus, we find 



C(6 1 )a + S(9 1 )( = 0, 
0(9^(3 + 5(^)7 = 0, 
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where C and S are defined in (|1.12l) . This can be rewritten as : 
(4.54) C = Ta and 7 = T(3 , 

with T defined in (Ti~Tll . 

4.2.2 Hessian matrix of (3 

A computation associated with the relations (j4.54l) gives, for the first deriva- 
tives at Xq\ 

a' 



dr\\(3\\ = (a cos 9 + Csin^) = a— 



dj/3|| = (/?cos0 + 7 sin0) = /^, 

and for the second derivatives ar xq\ 

2 

d 2 ||/3|| = 5 2 + ^-a 2 + 2D cos 9-2E sin9, 

O 

2 

d 2 s\\P\\ = e 2 + ^-p 2 + 2F cos9 - 6H sin9 
2 

d 2 rs \\/3\\ =Se + a(3^ + C cos 9 - 2B sin 9. 

So, we get the following Taylor expansion : 

/ / 2 2 

\\(3\\ = 1 + a^r + p^-s + T —{5 2 + ^a 2 + 2D cos9 - 2E sm9) 

22 2 
+ S — (e 2 + ^/3 2 + 2F cos 6> - 6H sin 6>) + rs(5e + ^a{3 + C cos 9 - 2B si 

We now want a Taylor expansion of 9. We have : 

sin 9 = sin9 1 + (cos 9 x d r 9 + sin 0^0) + r 2 {^^d 2 ^ - ^^(^) 2 ^ 

+ s 2 (^^d 2 9 - S -^l(d s 9) 2 ^j + rs (cos 9 1 d rs 9 - sin 9 1 d r 9d s 9) . 
After computations using (14.531 ). we get : 



d 2 9 2 aa' sin9 
— -a-^-DosmO-EoccsO- — 

d 2 9 ^aa' _ . . __ . sin# 



sin# 



d 2 9 = -2a(5^- - C sin 9 - 2B cos 9 - ^e<5. 

b z cos 9 
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Finally, we can write the second derivatives of (3 



(4.55) \& = ^fa 2 + 5" + - WW 

(4-56) ^ 2 /3 = ^tff fP + + F C(£) - 3tf o S(0) 



(4.57) d 2 J = ^ n t,Cf <*P + ^oSe + C C(8) - 2B S(6) 



2 sr ^ 2S(6) 2 " 2cos# 

- 3 + rV' , , g(g) 

2S(6>) 2 P 2cos0 ( 

5 Construction of a quasimode 



In this section, we construct a quasimode and prove Theorem 11.51 More 
precisely, we proceed in three steps. First, we choose an appropriate scaling. 
Then, we make a Fourier transform and a translation as annonced at the 
beginning of Section 2. We expand the new operator in powers of h and look 
formally for a quasimode expressed as a series in powers of h. As we are 
just interested in a three terms asymptotics for the bottom of the spectrum, 
we will just keep the three first terms of the series. The construction of 
the quasimode (denoted by if)) will then deeply use the results of Section 
2.2. In particular, the solutions f ,w ,d s ug and d 2 u g of equations of the 
form h(6)v = f will appear in the expression of the quasimode. Finally, we 
will estimate a quantity expressed as \\(P^ — A(/i))x^|| (where x is a cutoff 
function) by using the exponential decay properties proved in Section 2.1 to 
control the error of the truncations and we will conclude with the spectral 
Theorem. 

5.1 Rescaling 

We first make the following rescaling : 

(5.58) t = h 1/2 t, s = h 1/2 §, r = h 1/2 f. 

We omit the tilde, we divide by h and, instead ofH , we consider : 



(5.59) H*>™ = < pp D t m h app D t + /* (D r + + f 2 h 2 (D s + A 

+ fUD r + X )(£>. + A h s ) + fUD s + A h s )(D r + X) 
+ hR 1 (r,s,t)D r + hR 2 {r,s,t)D s , 
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where 



A r = V 6l + hV 2 P r , 2 + hP r , 3 , 
X = h l ' 2 P s , + hP s , 3 , 



and with 



n 



app 



1 + h 1/2 tK^ - /iTr(^i) + hR, m h app = 1 - h x ' 2 tKf + hTr(^±) + hR 



= 1 + 2h 1/2 tK n - hG n - hR 
= 1 + 2h l ' 2 tK 22 - hG 22 - hR 



app - ■ - u 1 v o 

/* = 2h^ 2 K 12 t - hG 12 - hR 12 



5.2 Quasimode 

We begin first by taking a partial Fourier transform in r, denoted by T and 
then we perform the translation Uq : 

T 

s = s + 



sin 



rh,new 



to get the operator T-C"-' new which is unitarily equivalent to H 

H h,ne W = u^T^H^^TUe. 
So, to get 7-L h ' new 5 we just have to replace (omitting the tilde) in H h,new 

T 1 

s by s H and r by D T -D s . 

sin sm 6/ 

We formally expand H h > new in powers of h 1 ^ 2 and we write : 



n 



h,new 



with 



H = D 2 + D 2 + V? 



Fx = (2V e (at - (s) + 2<5tL> s ) D T 



2C 
sin # 



+ (-^(^ - 7*M + -Aa(V B D. + D s V e ) + —JD s )\t] - -\-V t 
sm 6 sm 6 smf — sm 6 
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where H\ is defined by : 

Hi = ~—Mt - (s)(D s V e + V B D S ) + 2V e (/3st - ]-s 2 ) + V V e t 2 
sin u Z 

- — a Vl + et(sD s + D s s) + it 2 D s + H K 
sm U 

and where H K is defined in (11.131) . The reason why we have boxed the terms 
involving r is that, after a partial scalar product in the variable (s,t), we 
will be reduced to operators in the variable r. The expression of H 2 will be 
discussed later. 

We look for a quasimode expressed as a formal series in powers of h 1 / 2 : 



+00 

= Y, hj/ \ 

j=0 



attached to an eigenvalue in the form : 

+00 

j=0 

Then, writting formally that u h satisfies 

n h ' new u h = xy i , 

we are led to the system : 

h° : HqUq = \ u , 
h 1/2 : Hiu + HqU\ = Xiu + X u±, 

h : H 2 u + + H u 2 = A 2 tio + ^iUi + Ao«2- 

In the next subsection, we construct a solution for this system. 
5.2.1 The coefficient of h° 

H (as an operator in the (s,t) variable) has been analyzed in Section 2. We 
take as Ao the lowest eigenvalue of H (which was denoted before by H{6)) 
and we choose 

(5.60) u (r,s,t) = (po(r)ue(s,t) , 

with ())q (of norm 1) to be determined later. 
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5.2.2 The coefficient of hh 

We have to compute : 

< HiUq, u e > . 
This quantity has the structure : 

+ ^rD T (j) + Ct(/>o + 2>r 2 o + <£</> , 

where : 

21 =< (2V g (at - (s) + 25W s )u e , u e >, 
2C 

23 = - < -r— ^VeUfl, M >, 
sm U 

€ =< (2((3t - r s)V e + -^-(V e D s + D s V e ) + -^-tD s )u e , u e >, 

sm 6 sm t/ 

7 

55 =< — r-9-7^W0,U0 >' 
sm" 5 6> 

<£ =< Hiu e ,u e > . 

Let us examine 21. We recall that (12.231) holds. Moreover, using Lemma l2.10l 
and (14.541) . we find : 

< 2Vg(at — C s ) u e, u e >— 0. 
Thus, we have proved that 21 = 0. 

Let us now consider 23 and 2). By Lemma [2T4l we have : 33 = 2) = 0. 
Finally, let us prove that <£ is also equal to 0. 
Again by (|4.54p and Lemma 12. 101 we find : 

< (2(f3t -~fs)V e u e ,u e >= 0. 

Then, as we have just proved above, we get : 

2e 

< ^—ztD s ue,ue >= 0. 
sm v 

Finally, we recall that (12.201) holds. We deduce that Ai = <£. Let us simplify 
its expression. We first observe that : 

< t 2 D s ug, uq >= 0. 

Then, we get : 

< (at - (s){D s V e + V D s )u , u e >= 0. 
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Indeed, an integration by parts gives : 

< (at - (s)d s (Vgu e ),ug >= - < (at - (s)d s u e , V e u e >, 
where we have used Lemma [2.41 Moreover, in the same way, we get : 

< (D S S + SD S )U 9 , Ug >= 0. 

Consequently, we have (cf. (jl.101) and relations (14.441) ) : 



E = Ai = 2 < Vg(f3st - j-)u g , Ug> +T]< t 2 VgUg, Ug > ^~ 

2 sin 8 



t\D s Ug\ 



t>0 



+ < H K Ug,Ug > . 

So, for ui, we take a solution of the equation 

(H - a(9))u 1 = \iu - HxUq. 
Gathering all the results of Section 2.2, we can write the expression of u± : 



. a 



ui = (—w + 2iSfo) \ D T (j)o 



+ 



sin 2 8 



d s u e 



tD 



t9o 



+ 







-Wo - 1 



2 sin sin 



TpO 



1 



2 sin d 6 



d s u e 



r 2 0o 



+ Ml, 



' S sin 9 
where u\ is defined by : 

(5.61) (H-a)u 1 = (\ 1 -H x )ue. 

Before analyzing the coefficient of h, let us state a lemma concerning u\ : 
Lemma 5.1 We have the identity : 

2 sin 9 / VgUiUgdsdt = I d s (HiUg)ugdsdt. 
Jt>o Jt>o 

Proof. 

We take the derivative of f!5.6ip with respect to s, we multiply by ug and 
integrate. 



At this step of the proof, we have not used yet the non-degeneracy condition 
(11.141) . In fact, this construction is enough for the proof of Theorem 11.41 
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The coefficient of h 

We are interested in A2 and u 2 and we are led to examine the compatibility 
condition : 

< H1U1 + H 2 u ,u e >= A 2 0o- 

We have not given yet the explicit expression of H 2 , but it is easy to observe 
that the last equation can be put into the form : 

{A 1 t 2 D 2 t + A 2 r 4 + A 3 r 3 D T + A 4 rD 3 r + A 5 tD 2 + A 6 r 2 D T + A 7 r 3 + A 8 D 3 T 
+ AgD 2 + A 10 tD t + A n r 2 + A 12 D T + A 13 t + A lA )(j> Q = A 2 o . 

Moreover, the coefficients of K and G\ (introduced in Section 13.21) just play 
a priori a role in the coefficients A* for % 6 {9, • • • 14} (for degree and scaling 
reasons). In the next paragraph, let us be more accurate on that point. 

"Reduction" to the case of the flat metrics 

We prove that the coefficients of t 2 , D 2 and tD t do not depend on Gi and K 
and that we can do, in some sense, exactly as if the metrics were flat. Before 
starting the analysis, let us recall that the relations (14.441) still hold so that 
the curvature play a role in £ and i]. 

Let us first consider the dependence of the coefficients on K. We begin by 
noticing that the one of D 2 has no such dependence. 

Then, we consider the one of r 2 . Collecting all the terms involving K in 
< HiU\,uq >, we find that they will appear only in the following coefficient 
(we use Lemma I5TT1 ) : 

- 7 o - / Hid s u e u B dsdt - - 7 o - / d s (Hiu e )u e dsdt. 
2sm 3 9j t>0 2sm 3 9j t>0 

Thus, we can compute the commutator H\d s — d s H\ modulo the terms which 
do not involve the curvature ; we get : 

— \-iKut 2 sin 9 cos 9 - AK u t sin 9V 9 - AK l2 t sin 6D S ). 
2 sin- 3 9 

Comsequently, the coefficient is : 

< - . 2n (Knt 2 cos 9 -4:K n tVo)uo, u e > . 
2 sin 9 

Then, we consider the coefficient of r 2 involving the curvature in H 2 : 

7^11 cos 9t 2 2jKutV e 
2 sin 2 9 sin 2 9 
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In conclusion, the coefficient of r 2 does not depend on K. 
We now analyze the coefficient of tD t . Modulo terms not depending on the 
curvature, we see that there is only one term in Hi which has to play a role 
in tD t and looking at the expression of the corresponding coefficient, we are 
reduced to the same computations as for r 2 . 

Let us now consider the dependence of the coefficients with respect to G\ 
(see Remark l3~2l) . First, let us consider the term in r 2 and let us look for the 
terms involving G\ : 

-G& [ V?uj-G&[ \D s u e \ 2 
Jt>o Jt>o 

and : 

(5.62) ~t~~2~2) \QH [ sVgu 2 + IF [ tV e u 2 A . 

sin 6 V Jt>o Jt>o ) 

A priori, in (15.62ft . there should be the term in F, but the coefficient is 
essentially : < tD s ue,ug >= 0. Recalling the relations (14.461) . (15.621) gives 
the contribution : 

-- -Te( g u + G i) sind I sV e u l + ~~ • 2~o cos flG^ / tV e u 2 e . 
an 6 J t>0 sin 6 J t>0 

Then, we recall Remark 12.271 and the coefficients of Gf t and Gf 2 cancel. 
Secondly, we look at the terms in D 2 . and tD t and we observe that the terms 

2 2 

G r n and G 22 on one hand and G\\ and G r 2 s 2 on the other hand cancel for the 
same reason as for G{ ± and G 22 after having used the expressions (l4.4Tf) and 
(14481) . 

Finally, the terms in R (see Remark 13.21) playing at the scale h and with 
degree less than 1 with respect to r and s will just provide, after the Fourier 
transform and the translation, coefficients of D T and r and, as we will see, 
these terms are not important to prove Theorem 11.51 

Now we give an explicit expression of H 2 (modulo terms which do not change 
the coefficients of D 2 , tD t and r 2 , but only those of r, D T and the constant) : 

H 2 = P r ,2 (d t - —^——-D s , s + t] + P S)2 (d t - -^.D s , s +-?—,t) 
\ sine/ sine/ / \ sine/ sine/ / 

+ P r)3 ( D T - ~^D S , S+^—,t)V()+ VeP r ,3 ( D T - ~^D S , S + t] 

\ sine/ sine/ J \ sine/ sine/ / 

+ F s , 3 \D T - ~^D S , 8 +-?—,t)Ve + V e P s ,3 [D T - ~^D S , a +-L-,t) 
\ sm c/ smw / \ sin U sm J 
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where the Py are given in (I4.49p (and in which we put the index at all 
coefficients) . The guess is that we have the following cancellations : 



Then, we would be reduced to an equation of the form : 

(A 9 D 2 T + A 10 tD t + A n r 2 + A 12 D T + A 13 r + A 14 )</> = A 2 0o- 

Finally, writing formally that D T = r and r = s, we would hope to get the 

expression of the Hessian matrix. 

Let us now prove the cancellations (15 .631) ! 

Coefficient of -D 3 , rD 3 ., 

Such terms do not appear. Indeed, with the choice of gauge in f!4.49p . there is 
no r 2 in the homogeneous polynoms of order 2 and no r 3 in the homogeneous 
polynoms of order 3. 

Coefficient of r 4 

We can observe that the coefficient of r 4 in H 2 is : 



(5.63) 



Ai = Wi G {!,-•• 



8}. 



4sin 4 #' 



moreover the coefficient of r 4 



in < H\Ui, uq > is : 



2 sin' 



< Vgd s Ug, Uq > . 



But, taking the scalar product of ( 12.201 ) by uq, we obtain 



< V s d s ug,u e > 



sin 9 



2 



Coefficient of r 2 D 2 T 

The coefficient of t 2 D 2 in H 2 is : 



c 2 



sin 2 6> 



and in < H\U\,uo > ■ 



2( 2 
sin 3 9 



< V e d s u e , u e > . 
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Coefficient of r 3 

The coefficient of r 3 in H\Ui is : 

\ a < (St + Cs)d s u e , V e u e > + ; s 5 < (D s V g + V 9 D s )d s u e , u e > 



sm 6/ d sm & 2 sm 6 1 

7 e 

37; < Vgf ,Ug > . 



sin 3 

fn H 2 , the coefficient reads : 



< (St + u e > + 2 ^ < V^, « e > . 



5 sin 9 sin 6* 

Using Lemma [2771 we find that the term in disappear. 

Let us consider the term in 7^. We take the derivative of (12.201) with respect 

to s and make the scalar product with ug to find : 

< Vg8 2 s Ug, Ug >= and < (D s Vg, VgD s )d s Ug, Ug >= 0. 

We treat the term in with Lemma l2.12l and the wished cancellation follows. 
Coefficient of r 3 D T 

From < HiUi,Ug >, we find the coefficient : 

7C 



sin 4 

and from < H 2 ug,ug > : 



< Vgd s Ug,Ug > 

Co 



2sin#' 

Then, we conclude exactly as for the coefficient of r 2 D 2 . 
Coefficients of tD 2 and t 2 D t 

Their cancellation also comes from (12.201) . Lemma \2. 71 and Lemma \2. 121 We 
now study the expressions of the coefficients of the quadratic terms t 2 ,D 2 
and tD t . 

Coefficient of D 2 

The terms in D 2 coming from < H\U\, ug > are : 
a 2 

— < 2Vg(Cs + St)w ,Ug > -A5 2 < td s fo,Ug > 

'Zibet XT - / — v /~» % /* 2iab 

+ < %Vg(Cs + St)f , Ug > < td s W , Ug > 
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By Lemma 12.131 the imaginary terms cancel. Then, the terms in D 2 in 
< H 2 Uq, ug > are : 

a 2 

— <(Cs + Stfu e , u e > +S 2 < t 2 u e , u e > 
+ 2D < tVgUg, ug > +2E < sV e u e , u e > . 

So, we observe, with Lemma [2T6l and (]2.31l) . that the coefficient of S 2 is 

sine r 2 , , r 2 2 , , i r T , 2 , , c(o) 

/ tsugdsdt + / t Undsat = / tVgUgdsdt = -. 

cos 9 J t>0 J t>0 cos 9 J t>0 2 cos e 

Then, with Lemma 12.111 and again (12. 31ft (for the coefficients of D and E ), 
we have proved (see (14.551) ) that the coefficient of D 2 is : 

\d 2 $(x ). 

Coefficient of r 2 

From < HiUi, u g > and with Lemma IBTTl we get as coefficient of r 2 : 



< ( -^(^ " T*)Vfl + -^2-MD s + D s Ve) + -^tD s 
sin 9 sin 9 sin 9 



6 sin 2 sin 3 6> sin / 2 sir 

From < H 2 u ,ug >, we find the coefficient : 

< (fit — JsYug, Ug > ^2~— < (/^ ~ 77 s2 + 7^ 2 ) U 0i % > 



sin £/ sm # 

C 2 C 
+ —47: < D 2 ug, ug > +— -3- < (D,(/?t - 7s) + (/% - js)D s )u e , w e > 
sin 9 sin & 

+ 2 ,„ < tVgUe, M e > +6 , 2 - < sV^ue, w e > +2F < t-D s ^, « e > 
sm 9 sm (/ 

e 2 7 

+ . 2 /1 < ^ w 0> M e > < ( at - (s)D s Ug, Ug > . 

sin 9 sm 6 1 

We observe that 2F < tD s ug, u g >= 0. Then, let us gather some terms. We 
have first, with the same formulas as for the coefficient of D 2 : 



(3 1 ln , 2 a 3 + cr 2 a' 

>——W Q ,Ue > +— — < (Pt-'ys) Ug,Ug >= 2 . 2 , 

S sm sin sin 9S(9) 2 
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There are other terms in (3 2 (or /3y) : 
(5.64) 



T T 2 T T 2 

- -t-jt; < (^ st ~ o s )u>e,ue > + . 3 „ < [2V&(/3at - -s ),d s ]u e ,u e > . 
sin c/ Z Z sin C7 z 

But, we have : 

PWst - |s 2 ), <9 S ] = 2 sin 9(J3st - |s 2 ) - 2((3t - 7 s), 

so, with ( 12.1 Op . we find : 

< [2F e (/3st- |s 2 ),<9 s ]w e ,u e >=< 2sinfl(/3st- |s 2 )u e ,u e > . 

Thus, ( 15.641) cancels. 

Let us gather the terms in 777 : 

(5.65) - — 1 \- < t 2 u e , u e > +;r^W < t 2 [V e , d s ]u e , u e >= 0. 

2 sm 6 2 sm u 

As terms in e 2 , we have (exactly as in the paragraph concerning D 2 with S 2 ] 

e 2 2 4e 2 f C(0) 2 
< i > — r^td s f ,u e >= , 2 -er. 



sm (/ sm (/ 2 sm c/ cos # 

Let us now consider the terms in £ 2 : 

c 2 c 2 

(5.66) < D 2 u e , u e > < W + d s V e )d 2 u e , u e > . 

sin t/ 2 sm 6 1 

After an integration by parts, we can apply Lemma [2751 and find that ( 15.661) 
cancels. Then, we have : 



< Woug, u e > +3^- < sV g u e , u e >= (F C(9) - 3H S(9)) . 

sm 6 sm 6 sm t/ 



Let us now prove that all the other terms vanish. 
So, let us first notice (Lemma 12.131) : 

4ie . . 2ze/3 

< —^{pt ~ 7s)V e Jo, u e > + < , 2 td s w , u e >= 0. 
sin t/ 6 sin 6/ 

Then, we have, with Lemma \2. 71 : 

< —PhiVeds + d s V e )f , u e > + < -r^-tcPug, u e >= 0. 
sm 6 sm t) 
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Let us write the remaining terms (for the last one, see Lemma [5.11 ) : 

2iQ iBC 
< -r-rniPt ~ ls)V e d 2 s u e , u e > + < , » A V e d a + d s V 6 )w , u e > 
sin 9 b sin 9 

C 

+ -tVh < (D s (3t - js) + (3t - -ys)D s )u , u e > 
sin 9 



7 27 
+ -r-^: < (at - (s)D s u e , u e > + , < [d s , (at - (s)(d s V e + Vgd s ) ]u e , u e > 
sm 9 2 sm 9 



After two integrations by parts, it is easy to see that : 

-r^- < (£> a (/3t - js) + (/3t - 73)D 8 )ti 9 , u e >= 0. 
sm t/ 

Then, using Lemma 12.121 we get : 



< -r^TiiPt ~ 7s)V e d 2 s u e , u e > + < . 3 „ (Vfrfl, + d a V g )w , u e > 
sm t/ b sin 6/ 

C/3 



25 sin 3 9 

We have : 



< (2C + 4(Cs + St)d s )u e , M e >= 0. 



< (at - C,s)D s ug, u e >^ 



sin d 6* 2 sin J 9 

It is easy to see that : 

< [8 a , (at - (s)(d s V e + V g d 8 )]u g , u e >=< (at ~ (s)[d s , (d s V e + V g d 8 )]u e , u e > 
= —2 sin 9 < (at — (s)d s ug, ug >= —( sin 9. 

In conclusion, we have proved that the coefficient of r 2 is : 

1 



2 sin 2 9 



d 2 s (3(x ). 



Coefficient of tD t 

From < HiUi,Ug > we get the following coefficient of tD t : 



< {2V {at-Cs) + 2StD 8 ){-^w o -i-4^^u e + 2i^-f o ),ue > 

b sm 9 2 sin 9 sm 9 



+ < (-^(^ " 7*)Vfl + -^r- Q (V d D s + D s ^) + J±-tD s )£w + 2i6f ),u e > 
sm 9 sm sm 9 b 

c 



sin 2 # 



< [if 1,5s] we, w<? > 
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From < H 2 u ,ug >, we find the coefficient : 

< (fist - -S 2 + -)ug, Ug > + - 2 < (at - (s)D s Ug, Ug > 



sin6 ^ 2 2 y ' sin 2 
2 2e<5 2 

< (at — Cs){Pt — JSjUg, Ug > H — ; — - < t Ug, Ug > 



sin 9 sin i 

^-(4B / sVgu 2 e dstdt + 2C / tV e u 2 e dsdt) . 
a Jt>o Jt>o 



sin 6* 

We do not have to redo all the computations ; indeed, these are exactly the 
same kind of computations that we have met before. We just have to recall 
the relations ( 14.541) and observe that all is divided by sin 9 and we get as 
coefficient : 

-t— d rs (3(x ). 
sm 9 

In conclusion, O satisfies the equation : 

(5.67) Gp(D T , -^)(p + A 12 D T (j> + A 13 r0 o + A 14 <f) = A 2 o . 

sm 9 

A priori A±2, Ai 3 and An are complex numbers. For (ci,c 2 ,d) G C 3 , we 
introduce : 

(5.68) <5 Aci , C2 ,d = &p(D T + Cl ,T + c 2 ) + d. 

Then, &/3 tCl ,c 2 ,d is just some shifted harmonic oscillator. After some (com- 
plex) translation to eliminate the linear terms, it is easy to put Equation 
(15.671) into the form : 

>C1,C2, 

for some complex numbers c±, c 2 and d. Thus, we choose for <j) the normalized 
eigenfunction associated to an eigenvalue X 2 = 7 n (€5 i g) + d. At this step, we 
do not know that d is real. We just know that \ 2 — d is real. 

5.3 Upper bound 

The computations of the last section lead us to choose as a quasimode : 

i) = X (r, s, t)v h (h~ 1/2 r, h-V 2 s, hT^H), 
where x is a smooth cutoff function in a neighborhood of Xo and where : 

V h = TUgU h , 

with 



u h = u + h xl2 u x + hu 



2- 



39 



To get Theorem 11.41 we would just have to keep the two first terms (which 
were obtained only under Assumption (jl.8D ). but we do not explicit the proof 
which is essentially the same as the one of Theorem 11.51 (see below). 
We recall (14.521) and we observe that (we use that v h has an exponential 
decay to control the commutators of x with the derivatives by 0(h°°)) : 

(5.69) h 2 J2\\(\ r \ 3 + \s\ 3 + t 3 )\D k D j tlj\\\ <Ch 2 h 3 / 2 h-\ 

(5.70) h\\(\r\ 3 + \s\ 3 + t 3 )\tlj\\\ < Chh 3 ' 2 , 

(5.71) /i^lK^ + s^t 4 )!^!)! < C7i/i 4/2 /r 1/2 , 

j 

(5.72) ||(|r| 5 + |s| 5 + t 5 )H|| < Ch 5/2 . 
Thus, we have 

\\P^-H M ^\\ < Ch 5/2 . 

Then, by construction of u h and by controlling the remainders in the expan- 
sion of H h,new with the exponential decay, it follows : 

\\(H M - (X h + K0 2 + { ln {& p ) + d)h 2 M\ L2{mappdrdsdt) = 0{tf' 2 ). 

Then, with the spectral Theorem, we have : 

d(a(Pl),X h + X 1 h 3 / 2 + ( 7 n(©/3) + d)h 2 ) = 0(h 5 / 2 ). 

Thus, as a{Pji) C I , A G 1, Ai G 1, we get necessarily that -y n (&p)+d G R. 
In particular d is real and we have proved Theorem 11.51 
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